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Abstract 

An analysis of soft breaking under the constraint of modular invariance is given. 
The role of dilaton and moduli dependent front factors in achieving a modular in- 
variant V so ft is emphasized. Further, it is shown that in string models tan j3 is 
no longer a free parameter but is determined in terms of a s tring and the other 
soft parameters by the constraints of modular invariance and radiative electroweak 
symmetry breaking. The above framework is then used to analyze the neutralino 
relic density consistent with the WMAP data at self dual points in the Kahler and 
complex structure moduli. One finds that the combined set of constraints arising 
from modular invariant soft breaking, radiative electroweak symmetry breaking and 
WMAP lead to upper limits on sparticle masses for [i > 0. These limits are investi- 
gated for a class of models and found to lie within reach of the Tevatron and of the 
Large Hadron Collider (LHC). Further, an analysis of the neutralino-proton cross 
section shows that dark matter in these models can be explored in CDMS (Soudan), 
GENIUS and ZEPLIN. While the analysis is carried out within the general frame- 
work of heterotic strings, the possibility that the results of the analysis may extend 
to a broader class of string models, such as models based on intersecting D branes, 
is also discussed. 



1 Introduction 

Recently the Wilkinson Microwave Anisotropy Probe (WMAP)[1, 2] has measured the 
cosmological parameters to a high degree of accuracy leading to a very precise constraint 
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on the density of cold dark matter (CDM) in the universe. Thus the data suggests an 
QcDhih 2 in the range 



Here VLcdm = Pcdm/Pc where pcdm is the mass density of the CDM and p c is the 
critical mass density needed to close the universe, and h is the Hubble parameter in units 
of lOOKm/s.Mpc. The implications of this constraint were analyzed in Ref.[3] in the 
framework of mSUGRA and quite surprisingly it was found that relic density constraint 
of Eq. (1) allows for scalar masses to be in the domain of several TeV which may even lie 
beyond the reach of the LHC (For other analyses of WMAP data see Ref.[4, 5, 6]). This is 
so because radiative breaking of the electroweak symmetry occurs in part on the so called 
Hyperbolic Branch/Focus Point region (HB/FP)[7] which contains an inversion region[3] 
where typically the Higgs mixing parameter \i is relatively small, i.e., \p\ « m ,mi/ 2 
where m is the universal scalar mass, mi/ 2 is the universal gaugino mass at the GUT 
scale M G ~ 2x 10 16 GeV, and p is the Higgs mixing parameter which appears in the 
superpotential in the form fiHiH 2 . In the inversion region the lightest supersymmetric 
particle (LSP) is the neutralino with mass ~ \p\ and is degenerate with the next to the 
lightest neutralino and the light chargino. Consequently there is a lot of coannihilation 
in this region allowing for the satisfaction of the WMAP relic density constraints even for 
very large values of mo,mi/2. 

In this paper we show that the situation is drastically different in a class of string 
models because of the more constrained parameter space [8]. Additionally the constraint 
of radiative breaking of the electroweak symmetry is more stringent allowing for a deter- 
mination of tan/3. Our underlying procedure is to consider effective low energy theory 
below the string scale for a heterotic string model[9]. We assume this theory to pos- 
sess a T-duality invariance, specifically an SL(2,Z) modular invariance associated with 
large-small radius symmetry. Thus we assume that the effective scalar potential in four 
dimensions depends on the dilaton field S and on the Kahler moduli fields Tj (i=l,2,3) 
and later we will also extend our analysis by including the dependence on the complex 
structure moduli C/j. We require the scalar potential of the theory to be invariant under 
the modular SL(2, Z) T transformations given by 



VcDMh 2 = 0.1126 
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(aidi - biCi) — 1, (ai,bi,Ci,di e Z). 



(2) 



We follow the approach developed in Ref . [8] (For previous analyses on soft susy breaking 
in string theory under the constraints of modular invariance see Refs.[10, 11, 12, 13] and 
for its phenomenology see Refs.[14]). Thus following the supergravity approach adopted 
in Ref. [15] we assume the superpotential to be composed of a visible sector and a hidden 
sector so that W = W v + Wh where supersymmetry breaks in the hidden sector and is 
communicated to the visible sector where it leads to generation of soft breaking. The 
Kahler potential is in general given by 

d = -ln(S + S + ^lb S? S log{Ti + fi)) - ]T log(T % + T t ) + K u Q\Qj + HjjQjQj (3) 

i i 

where 5f s is the one loop Green-Schwarz correction to the Kahler potential[16] and Q 
are the matter fields consisting of leptons, quarks and the Higgs. For the analysis in this 
work we assume 8f s = 0. The resulting moduli potential have a form which arises in a 
class of no scale supergravity [ for a review see [17]]. For the visible sector we assume 
a general form W v = ft-uQiQj + ^ijkQiQjQk- The outline of the rest of the paper is 
as follows: In Sec. 2 we give a pedagogical description of modular properties of various 
quantities that appear in soft breaking. We emphasize the role of dilaton and moduli 
dependent front factors in achieving a modular invariant soft breaking and an explicit 
proof of modular invariance of V so ft is also given. In Sec. 3 we discuss radiative breaking 
of the electroweak symmetry in the context of string models and show that tan (3 is not a 
free parameter but a determined quantity in string models. In Sec. 4 we use the modular 
invariant soft breaking and the determination of tan (3 to compute the neutralino relic 
density in the range allowed by the WMAP data. The cases fx > and /i < are found 
to have significantly different behaviors. For the case \x > (which appears to be the 
preferred value of /i[18] from the — 2 data[19]) it is shown that the modular invariant 
soft breaking, the string determined value of tan (3 and the WMAP constraint combine to 
produce upper limits on sparticle masses which appear to lie within reach of the hadron 
colliders for a class of string models. Implications for the direct detection of dark matter 
are also discussed. Conclusions are given in Sec. 5. 
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2 Modular Invariant Soft Breaking 



In this section we give a pedagogic discussion of the modular properties of the soft param- 
eters and also show explicitly that the soft breaking scaler potential is modular invariant. 
We begin by defining modular weights. Suppose we have a function f(Ti, %) which trans- 
forms under modular transformations of Eq. (2) so that 

f(T t ,f t ) -> f(T>X) = (icT t + dT^-icT, + rfr/C^T.) (4) 

then it has modular weight (721,722). For example, from Eq. (2) one finds that (Tj + 
Tj)— >(T/ + T/)= (iCjTj + di)~ l (—iCiTi + dj) _1 (Tj + %) and thus (Tj + %) has modular 
weights (—1,-1). Now G defined by G = d + ln(WW^) is invariant under modular 
transformations which implies that W has modular weight (—1,0) and has modular 
weight (0,-1) under the modular transformation of Eq. (2). Thus WW^ has modular 
weight (—1, —1) and \W\ has modular weight (-1/2,-1/2). Defining 9 W so that 

W = \W\e i9w (5) 

we find that e l6w has modular weight (-1/2, 1/2). The Dedekind function 77(2^) (where 
r}(Ti) = g 1/24 n„(l - q n ), q = e 2niTi ) has modular weight (1/2,0) while G 2 (7i) defined by 

G. l( T, = M + * (6, 
oTi Ti + Ti 

has modular weight (2,0). More generally if we have a function Q(T;) with modular weight 
(t2q,0) then the covariant derivative defined by 

transforms as 

D Ti [lnQ(T t )} -> ^[/ng^')] = (icT} + d) 2 ,D Ti [ZnQ(T t )] (8) 

and thus has modular weight (2,0). 

In Table 1 we summarize modular weights of several quantities that appear in the 
analysis of soft breaking. We consider now the specific model of Ref.[8]. For W h we 
assume that the modular weight is carried by the Dedekind function r)(Ti) so that Wh = 
F(S,Ti)/l\i =1 ri(Ti) 2 , where F(S,Ti) is modular invariant. For the Kahler potential we 
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consider the specific form d = D(S, S) -X)jZn(Tj + Ti)+J2i a (Ti + Ti) n%a C^C a , where C a 
are the charged fields with modular weights n l a . We introduce the notation 7 S = (S + 
S)G,S/V3 = hs\e lds , 7t, = (Ti + f^G,Ti/y/Z = IttJ^. The condition that vacuum 
energy vanish then takes the form 

I7s| 2 + El7rj 2 = l (9) 
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Table 1: A list of modular weights under the modular transformations of Eq. (2) 
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The soft breaking potential for this case is given by 

Vsoft = m§ /2 E(l + 3 E <| 7Ti | 2 )ct Ca + (£ B%w§ + £ A% y w^ + H - c -) ( 10 ) 

o j=l a/3 a/37 

where the normalized fields c a and m 3 / 2 (= e _G / 2 ) are modular invariant and w^l = 

(3) 

H a j3C a Cp and w a p = Y aj3l C a C^C^. Following the standard procedure the soft breaking 
parameters A and B° may be expressed in the form 

D/2-ie w 

= -V3m 3/2 [|7s|e- ifls (l - (S + S)d s lnY aPl ) 
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+ E l7rje-^(l + < + ^ + n\ - (7J + T z )d T MY aPl - (T 4 + T.K^G^))] (11) 



i=i 

and 



5 Q/3 = -m 3 / 2 ^ [1 + V3|7s|e- i9s (l - (S + 5)9 s Zn^) 
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+v^E iTrJe"^ (1 + < + ^ - (Tj + %)d T ln^ - (T, + f i )n* a/3 G 2 (T i ))] (12) 



i=i 



The modular properties of A Q a p and are not so manifest as written in the above 
form. This is so because dr^lnY^) and ^(/n/ia^) have no well defined modular weights. 
Indeed they are mixtures of terms with different modular weights. Thus we may express 
them as follows 

d Ti (lnY a ^) = D Ti {lnY aPl ) + (1 + nj,^)(T i + T,)" 1 
d Ti (ln^) = D Ti (ln^) + (1 + n^pXTi + f^ 1 

where Dr^lnY^) and DT^lnHap) have modular weight (2, 0) while (Tj + Tj) _1 has modu- 
lar weight (1, 1). Thus we see that dr^lnY^) and (Zn// Q| g) are a linear combination of 
quantities with modular weight (2,0) and (1,1). Using Eq. (13) we may express Eqs. (11) 
and (12) in the form 

D/2-i0 w 3 

= ™3/ 2 t^v^E bill Pi + K^G^TA + D Ti lnY aPl ] 

v / i=i 

D/2-ie w 

-m 3/2 VShs\e- ies (l ~(S + S)d s lnY aPl ) (14) 
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and 

D/2~i9 w 3 

B% = m 3/2 -=- {-1 + v/3£ | 7 tJ(T, + T,)e-^ [n^G^T,) + D Ti ln^}} 

v J i= i 

e D/2-ie w 

-m 3/2 V3\7s\e- t9s (l -(S + S)d s ln^) (15) 

The modular properties of A and B° can be easily read off from above. Using Table 
1, we find that in Eq. (14) the front factor e D /2 ~ i9w / y/f has a modular (1,0) and the 
remaining factors have modular weight (0,0) giving A an overall modular weight (1,0) as 
desired. The same also holds for B° as can be read off from Eq.(15). While expressions 
of Eqs.(14) and (15) are useful in exhibiting the modular properties of A and B°, they 
(as well as Eqs. (11) and (12)) seem to imply at least superficially that A and B° have 
a significant dependence on the modular weights. This is not really the case as we now 
demonstrate. To this end we define modular invariant couplings A a/ j 7 and Ji a $ so that 

Y aM = X^v'^^m), = ^ = ^-^^(TO (16) 
where A a/ 3 7 and ji a p have modular weights (0,0). We note the following identity 

d T JnY a/3l = d Ti lnK Pl - (1 + < + nj, + n 7 )G 2 (T,) + 1 + < + n £ + n i (17) 

which can be gotten by using the relation d Ti ln\ a/ 3y = D Ti (lnY a/3l ) + (1 + li^^G^Tj) 
and Eq. (13). Using Eq. (17) in Eq. (11) we may write A^p as follows 

D/2~i8 w 3 

A° af3l = m 3/2 + f t )e~^ [-G 2 (T t ) + d Ti ln\ aPl \ 

VJ j=i 

D/2-ie w 

-m 3/2 V3\ ls \e- lds (l - (S + S)d s lnY aPl ) (18) 

and a similar analysis gives 

D/2~id w 3 

B% = m 3/2 + v^E \1tM + Tje-^i [-G 2 (T t ) + d Ti lnjl a p\} 

VJ i=i 

D/2~ie w 

-m 3/2 V3\is\e- lds (l -(S + S)d s ln^) (19) 

In Eqs. (18) and (19) one finds that the modular weights have disappeared due to the 
cancellation arising out of the identity Eq. (17). Further, using Table 1 we find that the 
modular weight of A is manifestly (1,0) and similarly the modular weight of B° is (1,0) 
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as necessary to achieve a modular invariant V so ft- The importance of the front factors 
e D/2~t9 w j^j-j wag g rg ^. em phasized in the work of Ref. [8]. This factor is often suppressed 
or omitted in string based analyses of soft breaking. However, since it has a non vanishing 
modular weight, the modular invariance of V so ft cannot be maintained without it. Further, 
the front factor e D l 2 is also significant in another context as will be discussed in Sec. 3. 

The analysis given above is more general than of Ref. ([8]). We can limit to that re- 
sult by making a further dynamical assumption. Thus while A a/ g 7 and j2 a p have modular 
weights (0, 0) they could still have modular dependence through dependence on modular 
invariants. A significant simplification occurs if we assume that this dependence is iden- 
tical to that of the function F(S,Ti). Specifically we assume that c^A^ = dT i F{S,T i ) 
and dTifta/3 = dT i F(S,T i ). Under this assumption and also assuming that the potential is 
minimized at the self dual points T; = (l,e l7r / 6 ) one finds[8] 

A Q aPl = m 3/2 e D ^(3j2\ lT f - V3\ ls \(l - (S + S)d s lnY a ^) e -^) / fl (20) 

i=l 

B% = m 3/2 e D/2 -^(-l + 3E I7tJ 2 - ^7*1(1 - (S + S)d s ln^)e'^) / fl (21) 

i=i 

where f a arises due to the possibility of several degenerate vacua so that {f a } = 8, 4\/3, 
6, 3\/3. We can easily extend the above analysis to include the complex structure moduli 
Ui (i=l,2,3) which have SL(2, Z) modular transformations of their own under which the 
scalar potential is invariant. Including the Ui moduli the Kahler potential is of the form 

3 3 

d = -ln(S + S)-J2 log{Ti + - £ log(U t + Ui) (22) 
i=i i=i 

It is convenient to relabel the Ui moduli so that T 3+ i = U (i=l,2,3). Then assuming that 
the potential is minimized again at the self dual points for all Tj, Eqs. (20) and (21) 
are valid with i summed from 1 — 6 where f a assume the set of values 2"3 3 -f (n=0,..,6). 
Further, in Eq. (9) the sum over i runs from 1 — 6. Gauginos also acquire masses after 
spontaneous breaking of supersymmetry so that M a = \FLef a ~ l e~ G l 2 f aa (G~ l )lG b , where 
f a is the gauge kinetic energy function for the case of a product gauge group G = Y\Q a and 
has the expansion f a = k a S+E ^[C(g a )-EiT(R^ l )(l+2n^ l )-2k a S GS ]log(r ] (Ti))+. . . , 
where k a is the Kac-Moody level for Q a . In our analysis we limit ourselves to Kac-Moody 
level 1 and to the universal case so that 

m 1/2 = V3m 3/2 \ ls \e- ida (23) 
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3 Determination of tan/3 from Modular Invariance 
and EWSB Constraints 



Radiative electroweak symmetry breaking (EWSB) constraints are more stringent in string 
theory than in supergravity. Radiative symmetry breaking produces two constraints aris- 
ing from the minimization with respect to the vacuum expectation values of the two 
Higgs fields v\ =< Hf > and v 2 =< H% > of the minimal supersymmetric standard 
model (MSSM), i.e., dV/dv® = 0. Typically in supergravity one of these is used to de- 
termine fi and the other one to eliminate the parameter B in terms of tan/5. Now in 
string theory one has in principle a determination of the parameter \i and the fact that it 
should equal the fi determined via radiative breaking of the electroweak symmetry would 
be highly constraining. In practice there is no realistic string determination of /i (see, 
however, Ref.[20]), and thus pending such a determination we will continue to use the ra- 
diative symmetry breaking equation for the computation of the \i parameter. The second 
symmetry breaking constraint which eliminates B in favor of tan (3 (see, e.g., Ref.[21]) 
as a parameter in supergravity must be treated differently in string theory since one has 
a determination of B here. A remarkable aspect of the B determination is that it has 
a front factor of e D l 2 which is related to the string constant via the relation e~ D = $ 2 

9 string 

where g s trin g = hdi and fcj is the Kac-Moody level of the gauge group Qi and gi is the 
corresponding gauge coupling constant. Thus the second radiative symmetry breaking 
equation can be thought of as a determination of the string constant in the terms of pa- 
rameters at the electroweak scale. Alternately we can think of this second equation as a 
determination of tan /3 so that 

tan/5 = ^ p Mz+ _ m ^ )U (| - 1 + 3£ N 2 - V3| 7 s|(l - (S + 5)9^) I)" 1 (24) 
V^^m 3 / 2 r B a strin g { 

where ?b = ?"b(1 + ez) is the renormalization group coefficient that relates B at the 
electroweak scale to B at the unification scale so that B = r B B and ez = Ml/{ra 2 Hi + 
/i 2 — |M| + (m 2 Hi + /i 2 + |M|) cot 2 /?). In this analysis we adopt the procedure of using 
Eq. (24) to determine tan (3. 

However, care should be taken in implementing Eq.(24). Since Eq.(24) arises from 
minimization of Vmggs at a higher scale Q, i.e., Q ~ rrif or Q ~ (highest mass of the 
spectrum) /2 [7], where the correction to /x 2 from the one loop correction of the effective 
potential is small, one has to consider the value of m 2 Hv rn 2 H2 , and ji 2 at this scale Q while 
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using Eq.(24). This is important since mf ll and m 2 H2 can vary a lot between Q ~ Mz and 
Q ~ mi etc. In this analysis considering one loop correction to the effective potential we 
have adopted an iterative procedure to determine tan (3. This is done by starting with a 
guess value of tan f3 then using it to compute B (via B(Q)) from the radiative electroweak 
symmetry breaking equation and then comparing with the string determined value of Bq. 
Such iterations produce rapid convergence of the Bq value as obtained from the radiative 
electroweak symmetry breaking toward its string value. Thus tan/3 is determined for 
given values of the soft parameters, input moduli etc., i.e., the choice of the self dual 
point and a strin g. 

4 Analysis of Dark Matter and Sparticle Masses at 
Self Dual Points with WMAP Constraints 

We discuss now the numerical results of the analysis. For simplicity we will set n l a = 
0, assume no dependence on the moduli Ui (the dependence on the C/j moduli will be 
considered later) and further assume that one is at the self dual point given by f a = 8. 
In the analysis we set all the CP phases to zero. In Fig. 1(a), for /i > 0, we give a plot 
of the self consistent determination of tan/5 and exhibit it as contours in the 7 S — m 3 / 2 
plane for values of m 3 / 2 ranging up to 2 TeV. One finds a steady decrease of tan/3 as 
7 S is increased for a given m 3 / 2 . Contours of A and ji are also shown. Regarding /i, 
for a given 7713/2 (i.e., for a given value of the universal scalar mass at Mq) /i increases 
with increasing 7 S because the universal gaugino mass at M G increases with increasing 75. 
Further, ji also increases with increasing 7773/2 for fixed 7 S since both the universal scalar 
mass and the universal gaugino mass at Mq increase with increasing 7773/2. In the analysis 
we also impose the flavor changing neutral current (FCNC) constraint from the process 
b -> s+7[22, 23, 24, 25] for which we take the range 2 x 1(T 4 < Br(b -> s+7) < 4.5 x 1(T 4 . 
In Fig. 1(a) the FCNC constraint is shown as a dot-dashed line below which the region is 
disallowed. Remarkably one finds that the WMAP relic density limits (0.094 < ^oh 2 < 
0.129) are satisfied within this very constrained system and most of the allowed parameter 
space is in the dilaton dominated region with j s > 0.6. The gray region-I in Fig. 1(a) 
refers to the the non-perturbative zone where the third generation of Yukawa couplings 
no longer stay within the perturbative limits because of very large tan (3 values. The gray 
region II results from the absence of radiative electroweak symmetry breaking (REWSB) 
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condition or smaller m-± below the experimental lower limit. In Fig. 1(b) we exhibit the 
satisfaction of the WMAP relic density constraints in the 7 S — tan (3 plane for fi > 0. 
The contours of m 3 / 2 are exhibited and so is the FCNC constraint from 5 ^ s + 7 (the 
disallowed region is in the left of the dashed line). From Figs. 1(a) and 1(b) one finds that 
consistent with WMAP constraint increases with increasing tan/5. However, even 
with tan/5 as large as 60 one finds that m 3 / 2 is no larger than 500 GeV in order that one 
satisfy simultaneously the WMAP constraints and the requirement of tan /5 becoming not 
too large so that all the Yukawa couplings stay within the perturbative domain. 

A plot of the mass spectrum as a function of 7 S with determined tan/5 but without 
the imposition of the relic density constraint is given in Fig. 2(a). A similar analysis as a 
function of m 3 / 2 is given in Fig. 2(b) but again without the imposition of the relic density 
constraints. Inclusion of the relic density constraint limits m 3 / 2 to lie lower than about 
500 GeV and the corresponding upper limits on sparticle masses can be read off from 
Fig. 2(b). In this case the upper limits on sparticle masses all lie roughly below 2 TeV 
and all are essentially accessible at the LHC. We note that the gluino is the highest mass 
particle over most of the allowed parameter space of the model. Inclusion of the FCNC 
constraint drastically lowers the upper limits so that the gluino mass lies below 500 GeV 
and masses of the other particles are even lower. Thus some of the particles may also 
be accessible at the Tevatron. Another interesting phenomenon is that over most of the 
allowed parameter space the neutralino is the lowest mass supersymmetric particle (LSP) 
and thus with R parity a candidate for cold dark matter. 

In Fig. 3(a) we give a plot of the neutralino-proton scalar cross section (7 y^p CIS 8b function 
of the neutralino mass for [i > 0. The region with black circles satisfies the WMAP 
constraint. Also exhibited are the sensitivities that will be reached by experiment, i.e., 
CDMS (Soudan) [26], GENIUS [27] and ZEPLIN[28, 29] (shown as two broad bands). In 
Fig. 3(a) we have included the limit from the first set of data from CDMS(Soudan) as 
announced recently [30] as well as the latest contour from EDELWEISS [31]. The dotted 
enclosed region is the contour from DAM A [32]. The analysis shows that in this class of 
string models dark matter falls within the sensitivity that will be achievable at the CDMS 
(Soudan), GENIUS and ZEPLIN. A similar analysis for the spin dependent neutralino- 
proton cross section a XP as a function of the neutralino mass is given in Fig. 3(b) for 
/i > 0. The region with black circles satisfies the WMAP constraint. 

The pi < case is analyzed in Figs. (4)-(6). In Fig. (4) contours of tan/5, and /i are 
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shown in the 7 S — plane. The constraint from b — > s + 7 is shown as a dot-dashed 
line below which the region is disallowed. The WMAP satisfied relic density region is 
shown as a small shaded area in black. The disallowed gray region I and III is discarded 
typically due to the absence of a tan/3 solution consistent with GUT scale inputs. The 
region II typically refers to the absence of REWSB, or having smaller than experimental 
lower limits of m~±. The region IV is a no solution zone like I and III, but its location 
and extent depends on the sensitivity of the minimization scale for REWSB. The region 
V is excluded because of tachyonic f±. An important point to note for /j, < case is that 
a given tan/5 contour has two separate regions, corresponding to smaller and larger 7^ 
values for a given m 3 / 2 . Furthermore, tan/3 is severely constrained to lie below 20 for 
/x < 0. Unlike the \i > case, the /x < branch has a larger region where |/x| is very 
small, close to its value constrained by the lighter chargino mass lower limit. Further, m 3 / 2 
here can acquire values much larger than in the /x > case consistent with the WMAP 
constraint and the REWSB appears to be occurring on the hyperbolic branch. Fig. 4(b) 
shows contours of m 3 / 2 in 7 S — tan/3 plane for /x < 0. The lightly shaded areas (in cyan) 
satisfy the 6-^5 + 7 limits. In contrast, areas marked with black dots do not satisfy the 
b — > s + 7 bounds. The WMAP allowed region is shown in small filled circles in black. 
Clearly, the WMAP constraint along with 6-^5 + 7 restricts tan (3 to be within 12 to 15 
. A plot of the mass spectrum as a function of 7 S is given in Fig. 5(a). A similar analysis 
as a function of m 3 / 2 is given in Fig. 5(b). Fig. 6(a) shows a plot of the neutralino-proton 
scalar cross section a XP vs neutralino mass for /x < 0. Here the detection cross sections 
are smaller than the case for /x > 0. However, the values of spin-dependent cross section 
of Fig. 6(b) for /x < are much higher than the case when /x > 0. 

We discuss now the dependence of the analysis on the choice of the self dual point. 
In Figs. (7) and (8) we give the analysis for the case when f a = 3\/3. A comparison of 
Fig. (7) with Fig. (1(a)) and with Fig. (3(a)) for the /x > case shows that the analysis at 
the self dual point f a = 3\/3 follows a similar pattern as the analysis at the self dual point 
f a = 8. The main difference is the lowering of the maximum in ms by about a 100 GeV 
without imposition of the 6 — * s + 7 constraint. However, when 5 ^ s + 7 constraint is 
imposed there is no significant change with respect to the result corresponding to f a = 8. 
For the /x < case a comparison of Figs. (8) with Fig. (4(a)) and with Fig. (6(a)) shows 
that the spin independent a XP cross section can dip lower for the case f a = 3\^3 compared 
to the case when f a = 8. Aside from that the features of Figs. (8) are very similar to 
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those of Fig. (4(a)) and with Fig. (6(a)). 

We now analyze the effect of including the complex structure moduli Ui. In Fig. (9 (a)) 
we present a composite analysis for relic density and Br(b — > s + 7) in (7 S — m.3/2) plane for 
/i > 0. Here we include Ui moduli and we scan over the self dual points corresponding to 
all possible values of f a : 2™3 3 ~ ri//2 , n = 0, ...6. This generates maximally allowed regions 
since we are integrating over all allowed values of f a . The b — > s + 7 contour is shown as a 
dot-dashed line below which the region is maximally disallowed. The discarded region by 
b — > s + 7 constraint extends up to m 3 / 2 = 750 GeV. The WMAP satisfied relic density 
region is maximally shown as small shaded area in black and the region is not much 
different from the result of Fig. (1(a)) and Fig. (7(a)). Similarly, Fig.(9(b)) does not show 
much different spin independent LSP-proton cross sections in comparison to Fig. (3 (a)) 
and Fig. (7(b)). We do a similar analysis for relic density for /j, < in Fig. (10(a)). The 
result shows that the WMAP allowed region lies in the range 7 S = 0.1 — 0.2 which 
is more constraining towards smaller 7 S values in comparison with the similar cases of 
Fig. (4 (a)) and Fig. (8 (a)). The spin independent detection cross section result is shown 
in Fig. (10(b)). In this case most of the points consistent with WMAP constraint lie at 
values of m x o » 500 GeV and fall outside the displayed range. Thus the number of point 
consistent with the WMAP constraint lying below m x o = 500 GeV is very sparce as is 
evident from Fig. (10(b)) where a few scattered points shown in black circles satisfies the 
WMAP constraint. However, the range of detection cross sections for WMAP satisfying 
points is not different from those of Fig. (6 (a)) and Fig. (8(b)). Thus we conclude that the 
inclusion of the complex structure moduli Ui does not cause any significant alteration in 
sparticle mass limits or detection cross section than the results without them. Another 
interesting observation is that the case \x > shows that the region consistent with the 
WMAP constraint is the dilaton dominated region where |7 S | > 0.5 while the case fi < 
shows that the region consistent with the WMAP constraint is the moduli dominated 
region where |7 S | < 0.5. Inclusion of U moduli for /i < further restricts 7 S so that 
|7 S | < 0.2. It should be interesting to see if this phenomenon is a generic feature in a 
wider class of string based models. 

The analysis above shows that dark matter that results in string models with deter- 
mined tan (3 is much more stringently constrained than in models where tan f3 is taken as 
a free parameter (For a sample of recent analyses see Refs.[33, 34, 35, 36, 37] and for a 
recent review see Ref. [38] ). While the above analysis is carried out in the framework of 
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heterotic string models, similar considerations may apply to other classes of string mod- 
els. Thus recently there has been a great deal of activity in the intersecting brane models 
(see, e.g, Ref. [39, 40] and the references therein) and specifically semi realistic models 
with 3 generations and N=l supersymmetry have been constructed[41]. The pattern of 
supersymmetry breaking in a broad class of such intersecting D brane models has been 
analyzed in Ref. ([42, 43]) and some common features with the soft breaking in heterotic 
string models are shown to emerge. Further, one expects that the constraint of radiative 
breaking of the electroweak symmetry in the intersecting D brane models would be very 
similar to the one in the heterotic string models allowing once again a determination of 
tan p. Thus we also expect upper limits on sparticle masses to emerge from the WMAP 
relic density constraint for the intersecting D brane models at least for the ji positive case. 

5 Conclusion 

In conclusion, in this paper we have analyzed the implications of modular invariant soft 
breaking in a generic heterotic string scenario under the constraint of radiative breaking 
of the electroweak symmetry. We have emphasized the importance of dilaton and moduli 
dependent front factors which are shown to be essential in achieving a modular invariant 
V so f t . Several forms of the soft parameters A and B° are exhibited, including the forms 
which exhibit explicitly their modular weights. It is shown that in models of this type 
tan /3 is no longer an arbitrary parameter but a determined quantity. The constraints of 
modular invariance along with a determined tan /3 define the allowed parameter space very 
sharply. Quite interestingly one finds that this parameter space allows for the satisfaction 
of the accurate relic density constraints given by WMAP. Further, the WMAP constraint 
combined with the FCNC constraint puts upper limits on the sparticle masses for the 
case n > which are remarkably low for a class of models implying that essentially all 
of the sparticles would be accessible at the LHC. Quite remarkably some of the sparticle 
spectrum should also be accessible at the Tevatron. Further, an analysis of the neutralino 
proton cross section indicates that this cross section should be accessible at the future 
dark matter experiments. The analysis presented here reveals that for the class of models 
considered the region of the parameter space consistent with the WMAP constraints is 
dilaton dominated for ji > and moduli dominated for ji < 0. It should be interesting 
to see if this result holds for a broader class of string based models where tan ft is again 
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fixed by the constraints of duality and the radiative breaking of the electroweak symme- 
try. Finally, we emphasize that the analysis presented here is the first work where the 
constraint of a determined tan f3 arising from the dual constraints of modular invariance 
and radiative breaking of the electroweak symmetry is utilized for the analysis of sparticle 
masses and dark matter. The predictions of this analysis on sparticle mass upper lim- 
its on the size of the neutralino-proton cross section have important implications for the 
discovery of sparticles at the Tevatron, at the LHC and for the discovery of dark matter 
via direct detection. Further, it would also be interesting to explore the implications of 
this predictive modular invariant scenario for other low energy phenomena, and for the 
indirect detection of dark matter including possible future experiments such as EUSO and 
OWL [44]. 
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(a) Contours of constant tan (3, Aq and \i are shown on the available parameter space for /i > 
in the (7^ — 17*3/2) plane. The b — ► s + 7 contour is shown as a dot-dashed line below which 
the region is disallowed. WMAP satisfied relic density region is shown as small shaded area in 
black. The gray region-1 refers to the discarded zone of very large tan /3 where Yukawa couplings 
are beyond the perturbative domain. The gray region II results from the absence of REWSB 
or smaller m~± below the experimental limit. 
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(b) The contours of constant TO3/2 in the 7 S — tan/3 plane for fi > are as labelled. The b — > s+7 
disallowed region is in the left of the blue dashed line. The WMAP allowed region is shown in 
narrow shaded area in black. 

Figure 1: 
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(a) Variation of sparticle masses with respect to 7 S when m 3 / 2 is fixed at 500 GeV for /i > 0. 
The WMAP constraint is not exhibited. 
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(b) Variation of sparticle masses with respect to TO3/2 when 7 S is fixed at 0.75 for fi > 0. The 
WMAP constraint is not exhibited. 

Figure 2: 
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(a) Scatter plot for spin independent LSP-proton cross section vs LSP mass for /i > when 7 S 
and are scanned. The region with black circles satisfies the WMAP constraint. Present 

limits (top three contours) and future accessibility regions are shown. 
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(b) Spin-dependent LSP-proton cross section vs LSP mass for /i > 0. The region with black 
circles satisfies the WMAP constraint. 
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(a) Contours of constant tan f3, and [i are shown on the available parameter space for /j, < 
in the (7 S — plane. The b — ► s + 7 contour is shown as a dot-dashed line below which 

the region is disallowed. WMAP satisfied relic density region is shown as small shaded area in 
black. The disallowed gray region I and III are typically discarded because of the absence of 
consistent GUT scale inputs. The region II typically refers to absence of REWSB or smaller 
than experimental lower limits of m-± . The region IV is a no solution zone like I and III, but 
its location and extent depends on the sensitivity of the minimization scale for REWSB. Region 
V is the tachyonic f\ zone. 
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(b) The contours of constant m.3/2 in the 7^ — tan /? plane for /1 < are as labelled. The shaded 
(cyan) regions satisfy the b — > s + 7 limits. The black dotted area does not satisfy b — ► s + 7 
limits. The WMAP allowed region is shown in filled circles in black. 

Figure 4: 
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(a) Variation of sparticle masses with respect to j s when m3/ 2 is fixed at 1000 GeV for \i < 0. 
The WMAP constraint is not exhibited. 
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(b) Variation of sparticle masses with respect to 7713/2 when 7 S is fixed at 0.35 for fi < 0. The 
WMAP constraint is not exhibited. 

Figure 5: 
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(a) Scatter plot for spin independent LSP-proton cross section vs LSP mass for /i < when 7 S 
and m3/2 are scanned. The region with black circles satisfies the WMAP constraint. Present 
limits (top three contours) and future accessibility regions are shown. 
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(b) Spin-dependent LSP-proton cross section vs LSP mass for /i < 0. The region with black 
circles satisfies the WMAP constraint. 
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(a) Analysis at the self dual point f a = 3-\/3 where contours of constant tan (3 are shown on the 
available parameter space for fi > in the (7=, — m.3/2) plane. The b — *■ s + 7 contour is shown 
as a dot-dashed line below which the region is disallowed. WMAP satisfied relic density region 
is shown as small shaded area in black. The gray region-I refers to the discarded zone of very 
large tan (3 where Yukawa couplings are beyond the perturbative domain. The gray region II 
results from the absence of REWSB or smaller m~± below the experimental limit. 
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(b) Analysis at the self dual point f a = 3\^3 where a scatter plot for spin independent LSP- 
proton cross section vs LSP mass for fi > is shown when 7 S and 7/13/2 are scanned. The region 
with black circles satisfies the WMAP constraint. Present limits (top three contours) and future 
accessibility regions are shown. 
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(a) Analysis at the self dual point f a = 3^/3 where contours of constant tan/3 are shown on 
the available parameter space for /i < in the (7 S — TO3/2) plane. The b — > s + 7 contour is 
shown as a dot-dashed line below which the region is disallowed. WMAP satisfied relic density 
region is shown as small shaded area in black. The disallowed regions I to V are same as those 
in Fig. (4(a)). 
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(b) Analysis at the self dual point f a = 3v3 where a scatter plot for spin independent LSP- 
proton cross section vs LSP mass for fi < is shown when 7 S and m 3 / 2 are scanned. The region 
with black circles satisfies the WMAP constraint. Present limits (top three contours) and future 
accessibility regions are shown. 
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(a) Composite analysis for relic density and Br(b — > s + 7) in (7 S — TO3/2) plane for /i > 
with the inclusion of Ui moduli at the self dual points corresponding to all possible values of 
f a : 2™3 3 -"/ 2 ,n = 0, ...6. The b -»• s + 7 contour is shown as a dot-dashed line below which 
the region is maximally disallowed. WMAP satisfied relic density region is maximally shown as 
small shaded area in black. The gray areas are similar to the discarded regions of Fig. (1(a)). 
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(b) Scatter plot for spin independent LSP-proton cross section vs LSP mass for fi > when 7 S 
and m 3 / 2 are scanned, for the composite analysis with (7; moduli corresponding to Fig. (9(a)). 
The region with black circles satisfies the WMAP constraint. 
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(a) Composite analysis for relic density and Br(b — > s + 7) in (-f s — 7713/2) plane for /j < 
with the inclusion of Ui moduli at the self dual points corresponding to all possible values of 
f a : 2™3 3- "/ 2 ,n = 0, ...6. The b — > s + 7 contour is shown as a dot-dashed line below which 
the region is maximally disallowed. WMAP satisfied relic density region is maximally shown as 
small shaded area in black. The gray areas are similar to the discarded regions of Fig. (4(a)). 
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(b) Scatter plot for spin independent LSP-proton cross section vs LSP mass for [i < when j s 
and m3/2 are scanned, for the composite analysis with JJ% moduli corresponding to Fig. (10(a)). 
A few scattered points shown in black circles satisfies the WMAP constraint. Most of such 
points corresponds to rh x a >> 500 GeV, thus falling outside the displayed limit. 
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